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Abstract. Several methods have been proposed in the literature to improve Free 
Electron Laser output by transforming the electron phase-space before entering the 
FEL interaction region. By utilising ‘beam by design’ with novel undulators and other 
beam changing elements, the operating capability of FELs may be further usefully 
extended. This paper introduces two new such methods to improve output from 
electron pulses with large energy spreads and the results of simulations of these methods 
in the ID limit are presented. Both methods predict orders of magnitude improvements 
to output radiation powers. 


PACS numbers: 41.60.Cr 

1. Introduction 

The Free Electron Laser (FEL) is an important scientific research tool that uses a 
relativistic electron beam to generate coherent radiation from the microwave through 
to the hard X-ray. At shorter wavelengths into the X-ray, this is unlocking many new 
areas of science in diverse fields such as: Warm-Dense matter studies |1]; short pulse 
protein diffraction [2] and medicine/surgery [3]. Current X-ray FELs [HIS] and those 
under construction [6], are unique laboratory sources of high power coherent X-rays. 
They are driven by electron beams generated from Radio-Frequency linear accelerators, 
which can be up to a few kilometres long. 

Many ideas are now being proposed to enhance and improve FEL output, towards 
shorter wavelengths, shorter output pulse durations, improved temporal coherence [7] 
and multi-colour operation [8]. These improvements extend the original high-gain 
FEL design where the electron beam from an accelerator is simply injected into a 
long undulator where the collective FEL interaction generates coherent output. The 
new methods rely upon manipulation of the electron beam in phase-space, using laser 
modulators and magnetic chicanes, either prior to injection into the FEL, or sequentially 
along the undulator as the FEL interaction progresses. 

Proposals also exist to reduce the overall lengths of FEL facilities by replacing 
the RF-linacs with plasma-wakefield accelerators laiini. These accelerators have large 
accelerating gradients about 10^ — 10'^ times larger than RF-linacs. However, the electron 
bunches generated so far are limited by a relatively large energy spread which inhibits 
any useful FEL interaction. As with the above proposed enhancements, methods 
that manipulate the electron beams have been proposed that may help mitigate the 
detrimental effects of energy spread. These include stretching the beam longitudinally 
before injection into the FEL to reduce the localised energy spread m, or transversely 
dispersing the electron beam to give a correlated transverse energy distribution and then 
matching this into a transverse gradient undulator [12] . 

Using a combination of modulators and chicanes, it is also possible to fourier- 
compose electron pulses of simple geometric shapes in longitudinal electron beam phase 
space e.g. rectangular, triangular, and sawtooth PI- Such waveform synthesis of 
the electron beam can also be utilised to generate phase-correlated harmonic beam 
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structures that can then perform analogous waveform synthesis of the coherent light 
emission from the beam structures. 

The electron beam parameters and manipulations described above can be very 
difficult, if not impossible, to model using conventional FEL simulation codes, which 
average the FEL interaction over a resonant radiation wavelength limiting both the 
radiation bandwidth that can be modelled and the range of electron energies, correlated 
or uncorrelated, within the beam. 

In this paper the nn-averaged FEL simnlation code PUFFIN na is used to simulate 
potentially usefnl electron beam undulator emission that would not be possible using 
conventional averaged FEL simulation codes. 

Firstly, a new method using electron beam phase-space manipulation is 
investigated, that may allow a FEL to operate with larger electron beam energy spreads 
which, for example, may assist the drive towards plasma-accelerator driven FELs. The 
method constructs a series of energy-chirped electron pulses (beamlets), each of different 
mean energy, vertically stacked in energy in phase space. The localised, or ‘slice’, energy 
spread of each beamlet is smaller than the original, unmodified beam from which the 
beamlets are constrncted. Previous work has used multiple beams generated individnally 
by a photocathode illnminated by mnltiple light pulses to generate different colonr pnlses 
from a FEL |T3]. Here, however, the beamlets are generated from a single electron pnlse. 

Secondly, an example is presented of what may be possible using fourier-synthesised 
electron beams [13]. This is the first simulation of the output from such waveforms 
in a FEL-type system. A fonrier-synthesised electron pnlse with a rectangnlar wave 
strnctnre in phase space is used to generate radiation in a series of nndulator-chicane 
modnles similar to those nsed in a mode-locked FEL amplifier [12]. The ‘discontinuous’ 
regions of the sqnare electron pulse form larger current regions that can emit significant 
coherent spontaneous emission (see e.g. m)- This coherent emission is periodically 
snperimposed using a sequence of undnlator-chicane modules and is shown to be able to 
generate significant output powers. This cannot strictly be called a FEL as little FEL 
interaction takes place. 

The methods simnlated here are clearly not to be considered as specific FEL design 
proposals, rather they are intended to demonstrate fnture possibilities and potential as 
electron beam generation advances beyond that of a simple linear beam model. 

2. Beamlets 

2.1. Beamlets - Description of Method 

In the Free Electron Laser (FEL), a relativistic electron beam of mean electron energy 
amplifies radiation in an nndulator of period A„ and rms magnetic field strength 
Bu. The resonant radiation wavelength amplified is given by = A^i(l -|- a^)/ 27 ^. The 
high-gain amplihcation process is characterised by the gain length Ig, where an initial 
radiation power Pq is amplihed exponentially as a fnnction of the distance through 
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the undulator as P{z) = Poexp{\/3z/lg) [IH]. With an electron beam energy of 7 ^, the 
gain length may be written, neglecting radiation diffraction and for no electron beam 
energy spread = 0, as: 

I = ^ ('!') 

^ Anp 2kuP^ 

where: = 27i/Xu, 
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is the is the FEL (or Pierce) parameter, au <x. Buku is the nndnlator parameter, oOp is 
the peak (non-relativistic) plasma freqnency of the beam, and Ipk is the peak cnrrent. 
For good amplihcation, the electron beam energy spread a-y mnst satisfy the ‘cold beam’ 
limit of: 


a, = ^ « 1. (3) 

Plr 

Optimal FEL gain is seen to occnr when Ipk is maximised and minimised. The 
method described below uses electron beam phase space manipulation to modify both 
of these parameters in an attempt to improve the FEL output potential of beams with 
large energy spreads (cTp >1). 

The method hrst generates a series of energy chirped beamlets stacked vertically 
in longitudinal phase space before they are injected into the FEL ampliher. As the 
FEL interaction occurs within the undulator further manipulation is required to ensure 
the radiation interaction with the chirped electron beamlets maintains a resonant 
interaction. 

In the hrst stage before injection into the FEL, the electron beam is passed through 
a modulating undulator and dispersive chicane, resulting in the beam phase space shown 
in hgure [T] This phase space is similar to the hrst modulator-chicane section used in the 
Echo Enhanced Harmonic Gain method [12]. The modulator-chicane sections perform 
the following consecutive transforms on the electron beam phase space coordinates: 

7 = 7o - Ay sin + 0^ (4) 

72 = 720 - 2D , (5) 

where the subscript 0 denotes the initial, untransformed coordinates, 72 = {ct — z)/lc is 
the coordinate in a window travelling at the speed of light scaled with respect to the 
cooperation length Ic = Xr/Airp of the FEL interaction. Ay is the energy modulation 
amplitude, n = Ai/A^ is the modulation period scaled with respect to the resonant 
wavelength and D = krpR^e is the scaled dispersive strength of the chicane. With this 
scaling, a resonant electron of energy y^ will fall behind a resonant radiation wavefront 
a distance Ic on propagating one gain length Ig through the undulator [20] . 
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It has been observed that in regimes where large dispersion is applied that the noise 
statistics of the macroparticles that simnlate the electrons in the dispersed beam can 
become incorrect. This occurs as the beam sampling in Z 2 is transformed into the 7 
dimension when rotated in phase space, and vice-versa. To ensure the correct noise is 
modelled, the functional form of the final electron beam phase space is used to initialize 
the beam before application of the noise algorithm isi and simulation using Puffin. 

A gaussian distribution for both dimensions of the initial beam phase space is 
assumed: 


/(^2,7) = 
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where: Zc is the electron pulse centre and are the standard deviations in 7 and Z 2 
respectively. 

By applying similar modulation and dispersive transforms to those outlined in 
the final beam distribution function obtained is: 
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Figure [T] plots the scaled longitudinal phase space distribution function of the 
electrons after the modulation-dispersive section and before injection into the FEL 
undulator using the scaled energy parameter pj = ( 7 ^- — 'jr)/p'yr with the following 
parameters: A 7 = 0.047^, D = 268.51, n = 68, cj) = 0, = 2p'jr (or dp = 2), 

■Jr = 1200, p = 1.6 X 10“^ and = 28.97. The modulation and dispersion of the beam 
is seen to create a stacked structure of energy chirped ‘beamlets’, slice sections of which 
are seen to have an energy spread which is reduced from the initial un-transformed 
beam with cXp = 2. Under certain conditions, each beamlet may then emit and amplify 
radiation independently of the other beamlets. The combined output from each of the 
beamlets may then give improved radiation output over the un-transformed beam. 

To illustrate how the method functions in the FEL undulator, a simplified version 
of the beamlet phase space is shown in figure]^ which consists of a series of chirped, zero 
energy spread, electron beamlets of different mean energies stacked in phase space. The 
chirp causes the radiation from one section of the chirped beam to drift out of resonance 
as it propagates into electrons which are resonant at a different wavelength. This 
impedes the FEL gain process. This effect may be successfully counteracted by using 
an appropriate undulator tapering to maintain the electron-radiation resonance m 
(These results have been reproduced using the simulation methods used here and are 
in very good agreement [23] •) Here, a different approach is demonstrated using a 
periodic series of undulator-chicane modules with multiple beamlets. The beamlets 
are periodically delayed by the chicanes so as to maintain a resonant interaction with 
the radiation generated by electrons of the same energy from the other beamlets. 
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Figure 1. The scaled longitudinal electron beam phase space distribution function 
given by equation Q f{z 2 ,p) (using p rather than 7 ) after transformation by a beam 
modulator and dispersive chicane. 


(Simulations using this method on the simple beamlet structure of hgure have been 
performed and presented elsewhere [23].) In the electron beam frame therefore, the 
radiation is passed from beamlet to beamlet so that it always interacts with electrons 
of a similar energy so maintaining a resonant interaction and giving an improved FEL 
interaction. This is achieved by making the slippage of a radiation wavefront through 
the electrons in each undulator-chicane module equal to the spatial separation of the 
beamlets. The enhanced slippage can also be expected to result in the generation of a 
series of modes in the radiation spectrum similar to that of [TB] which demonstrated that 
an undulator-chicane lattice will amplify side-band radiation modes that are separated 
by: 

Auj/ujr = dvrp/s, (8) 

where s is the slippage length in scaled units of Z 2 in one undulator-chicane module [16] . 

The FEL parameter p oc where Ipk is the electron pulse peak current, and is a 
measure of FEL efficiency. When considering individual beamlets a FEL parameter may 
also be defined for each beamlet: pb oc where Ib is localised (slice) current of the 
beamlet. (Note that as the beamlet energy is chirped, the mean pulse energy 7 ^, is used 
in the definition of pb-) Other beamlet parameters are also defined as pb = {'Jj — Jr)/Pblr 
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Figure 2. Scaled longitudinal phase space of the electrons for the simplified beamlet 
model. For a given energy the beamlets are separated spatially by Az 2 and for a given 
Z 2 the beamlets are separated by Ap = Aj/p^r- A chicane delay of the electrons 
corresponds to a positive shift in Z 2 - A series of chicanes slip the electrons forward in 
Z2 SO that they interact with the same resonant wavelength as emitted by the previous 
beamlet. 


and a beamlet scaled slice energy For a beamlet to lase independently its slice energy 
spread must then satisfy: 


^Pb 


a. 


lb 


Pb'7r 
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(9) 


(Note here, that the mean pulse energy 7 ^ is used in the definition of ap^ rather than 
a local ‘slice’ value 'jb. This can be considered a reasonable approximation for the 
inequality ([^, so long as 'ft does not differ significantly from 7 ^.) The beamlet slice 
energy spread cxp^ and instantaneous fractional FEL parameter pb/po, where po is the 
FEL parameter of the un-transformed beam, can be calculated and are shown in figures]^ 
and 1 ^ towards the higher energy and mid-sections of the electron pulse respectively. 

The energy spread condition for FEL lasing of equation ([^ may be used with the 
FEL radiation bandwidth saturation Au/ujr ~ 2p [21] to define the minimum energy 
separation A 7 of the beamlets so that the gain bandwidths of each beamlet do not 
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At the centre of the electron pnlse the beamlets split into pairs [25], i.e. two per half 
modnlation period, while for the electron pulse higher and lower energies, formed by the 
modulation extrema, the beamlet pairs merge into single beamlets as seen in figures 
andlH 

Both the energy spread condition ([^ and beamlet separation condition ( jloj ) are seen 
in £gure|^to be satisfied for the higher energy regions of the beamlets. (These conditions 
are also satisfied at the lower energy regions of the beamlets, but are not shown.) 
However, the condition placed on the beamlets’ energy separation (10) is not always 
satisfied at the pulse centre where the beamlets are formed in pairs, as seen in figure 
Hence, the undulator-chicane slippage length is set equal to the beamlet separation for 
the higher and lower energy regions of the pulse where the energy separation of the 
beamlets is approximately a constant. 

Results of a FEL interaction using an un-transformed (no beamlets) pulse with 
large energy spread cXp = 2 and of the transformed (beamlet) pulse are shown in figure]^ 
As expected, for the pulse without beamlets and the large energy spread, only small 
scaled peak powers of \A\^ ~ 10“^ are observed in the simulation. However, for the 
transformed pulse with beamlets that have smaller energy spread, (jp^ < 1 , and that are 
matched to the undulator-chicane modules, powers 2-3 orders of magnitude greater are 
observed. For the modulation period of 68 Ar used here [n = 68 ), matching was achieved 
using undulator modules of 20 periods and isochronous chicane slippages of 48A^. It is 
seen that the FEL lasing is greater for the lower energy beamlets of the pulse around 
Z 2 ~ 400. This preferential FEL interaction and amplification of the lower frequency is 
consistent with the scaling of the FEL parameter p oc 7 “^ which gives greater values 
and so strength of interaction, for lower beam energies. In the simulations here, the gain 
length of the higher to lower energy beamlets is up to ~50% larger. Evidence of the 
modal structure in the spectrum is also observed in the scaled power spectrum (inset), 
consistent with the undulator-chicane system which from ([^ gives a mode spacing of 
Aw = 0.0147. 


Signihcant bunching of the electrons in one of the lower energy beamlets, with a 
mean value of scaled energy < p >= (< 7 > —^r)lplr ~ —5, is also observed as shown 
in figure]^ Note from the lower plot for the spectrum that the electrons are bunched at 
a lower frequency w/wr ~ 0.85 than the mean resonant frequency of the electron pulse. 
This frequency shift from resonance is consistent with the lower mean energy of the 
electrons as Aw/w^ ~ 2p < p >= 0.16 and is in agreement with the radiation frequency 
spectrum of hgure Electron bunching is also observed in a higher energy beamlet of 
mean scaled energy < p 4, shown in figure Here, the bunching is seen to be at 
a less advanced stage, but can be expected to reach saturation on further propagation 
through the undulater-chicane lattice. 
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Figure 3. Top panel: Detail of the higher energy beamlet phase space distribution 
function of equation with a single beamlet delineated by white dashed lines. The 
values of the scaled energy spread (middle) and pb/po (bottom) were calculated 
for the single beamlet as a function of Z 2 - Towards the pulse head (z 2 < 105) the 
electron pulse is diffuse with a larger energy spread and smaller pb- Nearer the 
centre of the pulse (105 < Z 2 < 125), the scaled energy spread decreases as the local 
density, and pb increase. However, further towards the pulse centre Z 2 > 125 the energy 
spread increases further as the beamlet splits into two identifiably separate beamlets, 
while the value of pb tends towards a more constant value. The condition for lasing 
of the beamlet of ap^^ < 1 is seen to be satished within this the head of the pulse 
(and is also satisfied at the lower energy beamlets of the tail). The energy separation 
between beamlets is also seen to satisfy condition (10) so that each beamlet can lase 
independently. The energy separation between beamlets does not change significantly 
with Z 2 , as neither does the longitudinal separation of beamlet regions with the same 
energy. Towards the centre of the pulse however, the beamlet structures have a more 
complicated phase space structure. 


3. Fourier Synthesised Electron Beams 

Further types of phase-space transformation of an electron pulse prior to generating 
radiation have been proposed and called ‘beam-by-design’ [7]. An example is 
investigated here to demonstrate the potential of such beam transformation prior to 
injection into the FEL and the subsequent transformation in the FEL emission stage 
using a series of undulator-chicane modules. An electron pulse consisting of a series 
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Figure 4. As figure but around the centre of the electron pulse about the mean 
pulse energy. The beamlets are seen to ‘spilt’ into two separate beamlets. While 
the scaled energy spread requirement for lasing < 1, is satisfied, the beamlet 
energy separation condition (101 is only satisfied for a small region of beamlets about 
the pulse centre. The beamlets are therefore unlikely lase independently with non¬ 
overlapping bandwidths, so that the effective energy spread for the interaction is 
increased, decreasing the ability of achieving significant FEL lasing. 


rectangular shaped distributions in phase-space can be generated [13] and contains 
a periodic series of current ‘spikes’. These current spikes are a source of coherent 
spontaneous emission which may, through a series of periodic superpositions enabled 
by chicanes, generate signihcant radiation output from an undulator-chicane lattice. 
We note that other methods can generate similar beam structures, e.g. the E-SASE 
approach [2S], however the methods of [13] are used here to demonstrate the types of 
more exotic interaction that may be modelled using non-averaged simulation codes such 
as PUFFIN. 

3.1. The Model - Coherent Emission from Rectangular Electron Fulses 

A new approach to produce so-called ‘RF-function’ electron beams was introduced 
in [13] . An RF-function generator produces a series of repeated wave forms by combining 
sine-waves of different frequencies as in a Fourier series. In a similar fashion, an 
electron pulse can be created with a phase-space that consists of repeated ‘waveforms’ 
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Figure 5. A comparison of the scaled radiation temporal power and spectral power 
(insets) for an un-transformed electron pulse (panels a and b) and transformed pulse 
of beamlets (panels c and d), when propagated through an simple undulator and an 
undulator-chicane lattice respectively and interaction length of z « 30. The (red) 
box shows the position of the electron pulse relative to the radiation (the head of 
the pulse is to the left.) Note the different lengths of the electron pulses due to 
differing dispersive effects of the chicanes. The beamlets propagating through a simple 
undulator (panel c) is seen to give a small improvement to the output from the un¬ 
transformed beam through both an simple undulator and an undulator-chicane system 
(a and b respectively.) The improvement in output from the beamlets is increased 
significantly when they are propagated through the matched undulator-chicane lattice 
as shown in panel d. The undulator-chicane lattice amplifies side-band radiation modes 
generated by the undulator-chicane modules and are separated by Aw = 0.0147 as seen 
from the panel d inset and in agreement with the mode-spacing relation of ([^. For 
all results shown in this figure the radiation field has been filtered about the resonant 
frequency 0.5 < w/w^ < 1.5 to eliminate low frequency coherent spontaneous emission. 


by modulation the electron beam using a series of seeded undulator modulators using 
different seed wavelengths, amplitudes and phases. Following the notation of [13], here 
a rectangular beam shape in phase space using a triple modulator-chicane lattice is 
synthesized and subsequent radiation generation following injection into an undulator 
chicane-lattice is modelled using PUFFIN. 

While in [13] an inhnity long electron beam was assumed, here, a hnite electron 
pulse with an initial Gaussian distribution in both and 7 is assumed, as given by 
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Figure 6. Electron bunching in a lower energy beamlet at z « 30. The top panel plots 
the charge-weighted electron phase space distribution ; the middle plots the bunching 
parameter of the beamlet at the fundamental radiation frequency, and the bottom 
plots the bunching spectrum of the beamlet. 

equation (§. As detailed in the Appendix, the same Fourier synthesis as outlined 
in P! is applied using the beam modulation transforms given by equation Q and the 
energy dispersion transforms of equation ([^. 

In electron phase-space, the vertical segments of the rectangular waveform generate 
regions of enhanced current, albeit with a larger energy spread. Each period therefore 
contains two current ‘spikes’ which can generate significant coherent spontaneous 
emission when their width is of a similar scale to a resonant wavelength UTl. However, 
due to electron beam dispersion in the undulator, the sharpness of the current spikes 
reduce on propagation, resulting in diminishing coherent emission. This dispersion of 
the current spikes may be compensated for by the use of chicane systems with a negative 
dispersion to allow for more prolonged coherent emission. The design of chicane delay 
systems with negative dispersion have been previously designed and tested as part of 
an accelerator lattice [2Z] and are also necessary for generating the RF-function beam 
shapes mm- If the slippage per undulator-chicane module is also made equal to the 
current spike separation, then the radiation is propagated from spike to spike and, if 
correctly phased, can facilitate the constructive interference of the coherent emission 
from each current spike in each new undulator module. 
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Figure 7. Electron bunching in a higher energy beamlet at z « 30. The top panel plots 
the charge-weighted electron phase space distribution; the middle plots the bunching 
parameter of the beamlet at the fundamental radiation frequency; and the bottom 
plots the bunching spectrum of the beamlet. 


3.2. Results - Coherent Emission from Rectangular Electron Pulses 

The following simulations use the same electron pulse parameters as the previous 
section, i.e., the electron pulse’s large energy spread is prohibitive to FEL gain. 
The phase-space distribution of the electron beam for the rectangular waveform was 
constructed from the analysis of the Appendix for three undulator-chicane modules 
using the following parameters in [A'y,D]: [A 71 = Way, Di = nip 7 rA/ 37 r/( 2 A 7 i)]; 
[A 72 = A 7 i/ 4 ;T )2 = -3T>i]; [A 73 = A 72 / 16 ;T )3 = - 3 / 12/4 = 9/li/4], with ni, 2,3 = 20, 
01,2 = 0 and 03 = 71. 

The initial current prohle of the electron pulse contains a series of current spikes 
at half the modulation period corresponding to 10 resonant radiation wavelengths or 
10 X 47 rp 2 in units of Z 2 . On injection into an undulator, these spikes act as a periodic 
series of phase correlated coherent emitters which, for a relatively short interaction 
length of z ~ 1), generate a broad modal radiation spectrum as seen figure]^ However, 
it is seen that alternate current spikes have dispersed to leave a series of more prominent 
current spikes at twice the initial spacing of Az 2 ~ 4. This is reflected in the temporal 
separation of the larger radiation spikes separated by Az 2 ~ 4). This also agrees with to 
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Figure 8. The evolution of a rectangular electron beam in an undulator showing 
top-to-bottom, the scaled radiation power \A\'^ as a function of Z 2 , electron phase 
space ( 7 / 7 r,^ 2 ) with detail inset, the scaled electron current as a function of Z 2 and 
the logarithm of the scaled radiation power spectrum \A\'^ as a function of the scaled 
frequency w/w,.- The series on the left plot the output for a scaled distance through 
the undulator of z « 1 and on the right for z « 20.1 An electron pulse with an initially 
large energy spread has been transformed into an electron pulse that contains a number 
of rectangular waveforms (see second plot on the left). The electron pulse structure 
now contains a series of current spikes of spacing Az 2 ~ 4. When this electron pulse is 
passes through an undulator each current spike acts as a source of coherent spontaneous 
emission. The radiation spectrum (bottom panels) show a broad bandwith modal 
structure with modes separation Auj/uJr ~ 0.05. As the electron pulse propagates 
along the undulator, the rectangular waveforms will disperse, and increase the current 
spike widths and reduce current spike amplitudes. As the current spikes’ ’sharpness’ 
decrease the coherent radiation produced by the current spikes will decrease. Because 
of this no amplification is seen when passing such an electron pulse through a long 
undulator, as shown in the r.h.s. of this figure. 


the spectrum in which a series of modes are generated with separation, from equation ([^, 
of Au/ur ~ 0.05 about the resonant frequency. 

On propagating further through the interaction region to larger values of 2 ~ 20, 
the right hand panels of hgure show that the energy modulation of the rectangular 
electron beam causes the electron beam to disperse in the undulator degrading the 
visibility of the current spikes and so decreasing the coherent spontaneous emission 
generated. Clearly, these dispersive effects mean that there is no beneht in increasing 
the interaction length over that oi z = 207rp 1. 
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Figure 9. By using chicanes with a negative dispersion, the undulator dispersion 
of the rectangular sections of the electron beam can be partially compensated for as 
seen here for z « 10. In doing so the electron pulse can continue to emit coherent 
emission in each undulator module. Here, each undulator module has 5 periods and 
each chicane delays the electron pulse by approximately 5 resonant periods, to match 
the current spike separation. Note that there is small FEL interaction as evidence by 
electron microbunching (not shown). 


By using chicanes with a negative dispersion it is possible to partially compensate 
for the undulator dispersion and maintain a spiked current prohle that can continue to 
emit CSE over a larger number of modules. An example of this is shown in figure 
were the chicane dispersion is set equal the negative of the undulator dispersion, i.e. 
D = —I [23]. The total undulator-chicane slippage for the radiation was again set 
equal to the current-spike separation, s = 10 x Anp. For this case, undulator-chicane 
modules of 5 undulator periods and 5 chicane slippage periods were used. In this way, 
the CSE from successive undulator-chicane modules superimpose and constructively 
interfere increasing the radiation power emitted. 

However, the radiation helds from each undulator-chicane do not superimpose 
coherently and the radiation energy is seen (not shown) to scale approximately as the 
number of undulator-chicane modules - a phase-matched coherent superposition would 
give a radiation energy which scales as the square of the number of undulator-chicane 
modules. The reason for this non-coherent superposition is that the dispersion of the 
large energy modulated beam in the undulators cannot be perfectly compensated for 
by the negative dispersion in the chicanes. (Phase space dispersion of electrons in the 
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Figure 10. As figure but now using an optimised chicane which maintains the 
rectangular waveform electron pulse structure in phase space as it propagates through 
the undulator-chicane lattice. The rectangular electron waveform emits coherent 
radiation in each new undulator module which constructively interferes with the 
radiation in subsequent undulator modules. 


undulator is due to differences in the axial speed v^, while electron dispersion in the 
chicanes is due to differences in the electron energy, 7 .) This is observed from the slight 
‘bowing’ of the rectangular structure of the electrons in phase space in figure Two 
possible methods to improve this are to reduce the initial energy modulation of the 
rectangular wave (the results here are for a relatively large energy modulation) or to 
use a (hypothetical) optimised chicane design which has a non-linear dispersive strength 
as a function of 7 . Here the latter is used and the results shown in figure [T^ Now, 
the bowing of the rectangular structure of the electrons in phase space is seen to be 
removed and the power of the radiation increased. The coherent radiation from each 
undulator-chicane module is now phase matched and is superimposing coherently after 
each module. The radiation energy is now also observed to increase in proportion to 
the square of the number of modules. 

A comparison of a normal (untransformed beam) FEL amplifier with the methods 
of beamlets of the previous Section and that of the fourier synthesised rectangular beam 
of this section is given in figure which plots how the scaled energy E of the radiation 
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Figure 11. The total radiation field energy E{z) of the normal (untransformed, 
no chicanes) FEL, beamlet FEL, and rectangular beams with linear and optimised 
chicanes. The radiation is filtered around the resonant frequency, 0.5 < uj/uJj. < 1.5. 


pulses evolves with the interaction length z, where: 

/ +00 

\A{z,Z2)\dZ2. (11) 

■oo 

Before performing the integral in ( [IT| the field was first fourier bandpass filtered so that 
only contributions about resonance in the interval 0.5 < oj/ur < 1.5 are considered 
(this removes the signihcant low-frequency CSE content.) The introduction of the 
phase=space transform to generate electron beamlets is seen to increase the exponential 
growth rate over the normal FEL interaction by a factor of approximately two. While 
the rectangular electron beams are seen not to have an exponential gain, it is essentially 
a Coherent Spontaneous Emission process, the starting powers are much greater than 
the FEL processes which start from spontaneous shot noise. It should be noted that 
when the CSE simulations predict radiation powers that are a signihcant fraction of the 
electron beam energy, that the effects of photon recoil should be included in the model. 
These effects are not included in the classical simulations presented here. 

4. Conclusion 

This paper has sought to demonstrate what may be possible when electron beams are 
transformed to alter their properties before injection into an FEL-type system. It is 
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stressed that the methods demonstrated here are not proposals for any specihc design 
or operational wavelength. Rather, they are used to demonstrate possible research 
directions towards future light sources, some of which have already been envisaged [7]. 

Here, the focus was to generate signihcant radiation output from electron beams 
that have insufficient beam quality to lase under normal FEL operation. These methods 
may be developed further and made more specihc e.g. to the electron beams generated 
from plasma accelerator sources which, to date, tend to have relatively high energy 
spreads. Other possibilities, such as multiple frequency generation, ultra-short pulses, 
chirped pulses (possibly shorter wavelengths) and others, are potential research areas. 
One topic that is apparent, but has not been explored here, is the introduction of 
tapered undulators into the design process. For example, the introduction of tapered 
undulators, matched to compensate for the chirped beamlets of above, instead of using 
chicanes, can be expected to produce interesting radiation output. 

It is noted that the simulations presented here cannot be modelled effectively, or 
at all, using simulation codes that are used to successfully model the ‘normal’ types 
of FEL interactions. Unaveraged FEL codes, such as the PUFFIN code used here, are 
required. 
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Appendix 

The final distribution function of a triple-modulator-chicane scheme is given below. 


/(^2,7) = —— exp ( ^ 

2tt (Tz2 \ 2(ji^ 


[((7 + A73 sin(--(z2 + 2D-i{y - 'yr)hr) + bs)) 

2n3p 


+A72 Sin(-^((Z2 + 2L>3(7 - jr)/jr) + 2 £> 2((7 + ^73 Sin(--^(Z2 + 2£>3(7 - 7 r)/ 7 r) + bs)) - lr)/lr) + ^2)) 
2n2P 2n3p 

+ A 71 Sin(-^(((Z 2 + 2 D 3(7 - 7r-)/7r) + 2D2{i'-f + A73Sin(-^(z2 + - Jr)/lr) + h)) - lr)hr) 

Znip zn^p 

+ 2 Di{{{y + A73sin(-^(Z2 + 21)3(7 - 7 r-)/ 7 r) + da)) + A72 sin(-^((z 2 + 21)3(7 - lr)hr) 
2n3p 2n2P 

o\ 

+2£>2((7 + A73 sin(-- {Z2 + 2D3{y - -fr)/lr) + bs)) - 7r)/7r) + ^2)) - 7r)/7r) + bl)] “ 7r 

2 n 3 p 

exp (^[((^2 + 21 ) 3(7 - 7 r)/ 7 r) + 2 i) 2((7 + A 73 Sin(^^(z 2 + 21 ) 3(7 - 7 r)/ 7 r) + h)) - 7 r)l 7 r) 

+ 2 Di (((7 + A 73 sin(-^(Z 2 + 21)3(7 - 7 r-)/ 7 r) + da)) + A72 sin(-^((z 2 + 21)3(7 - 7 r)l 7 r) 
2n3p 2n2P 


+2D2((7 + A73 sin(-- {Z 2 + 21)3(7 - 7r)/7r) + da)) - 7r)/7r) + (j)2)) - 7r)/7r] “ 

2n3p 


The energy modulation parameters A 7 i 2 , 3 , modulation frequencies ^ 1 ^ 2,3 = 
kix^/kw and modulation phases 0 i, 2,3 are associated with first, second and third 
modulator sections respectively. Similarly £* 1 , 2,3 are the dispersion factors for chicane 
1,2 and 3. o '^^22 is the standard deviations in 7 and Z 2 . The resonant energy is defined 
as 7r =< 7 > | 2 =o and the electron pulse centre is given by 7c. 
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